library(sf)
package ‘sf’ was built under R version 3.6.2Linking to GEOS 3.7.2, GDAL 2.4.2, PROJ 5.2.0
library(tidyverse)
library(ggplot2)
loc_authority_with_extra_density_2019 <- read_csv("data/loc_authority_with_extra_density_2019.csv")
urban_rural_2019_with_ex_den <- read_csv("data/urban_rural_2019_with_ex_den.csv")
la <- st_read("https://opendata.arcgis.com/datasets/2698b6825e88404ab6272a5ff7f59ee6_0.geojson")
loc_authority_with_extra_density_2019
loc_authority_with_extra_density_2019 %>%
filter(feature_code == "S12000044")
loc_authority_with_extra_density_2019 <-
loc_authority_with_extra_density_2019 %>%
mutate(feature_code = recode(feature_code, S12000049 = "S12000046", S12000050 = "S12000044"
)
)
urban_rural_2019_with_ex_den <-
urban_rural_2019_with_ex_den %>%
mutate(feature_code = recode(feature_code, S12000049 = "S12000046", S12000050 = "S12000044"
)
)
la %>%
filter(lad18nm == "North Lanarkshire") #S12000046
Simple feature collection with 1 feature and 10 fields
geometry type: MULTIPOLYGON
dimension: XY
bbox: xmin: -4.19449 ymin: 55.73476 xmax: -3.711127 ymax: 56.03136
CRS: 4326
objectid lad18cd lad18nm lad18nmw bng_e bng_n long lat st_areashape
1 365 S12000044 North Lanarkshire 277984 665608 -3.9514 55.8681 472257663
st_lengthshape geometry
1 160705.9 MULTIPOLYGON (((-4.020121 5...
la_loc_authority_ex_den <- inner_join(loc_authority_with_extra_density_2019, la, by = c("feature_code" = "lad18cd"))
la_loc_authority_ex_den
la_urban_rural_2019 <- left_join(urban_rural_2019_with_ex_den, la, by = c("feature_code" = "lad18cd"))
la_urban_rural_2019